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We have investigated the two-photon nonlinearity at general cavity QED systems, which covers 
both weak and strong coupling regimes and includes radiative loss from the atom. The one- and 
two-photon propagators are obtained in analytic forms. By surveying both coupling regimes, we 
have revealed the conditions on the photonic wavepacket for yielding large nonlinearity depending 
on the cavity Q-value. We have also discussed the effect of radiative loss on the nonlinearity. 
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I. INTRODUCTION 

\q • The nonlinear optical response appears strongly when intense light fields are irradiated onto the material with large 
nonlinearities. By using an optical cavity, we can magnify the intensity of the input field inside the cavity, if the 
input field is resonant to the cavity mode pj. Thus, we can enhance the nonlinear optical response by putting the 
nonlinear material inside the cavity. This is realized experimentally by using two-level atoms as the nonlinear material 
@jES- This idea opened the possibility of obtaining large nonlinearity even by weak input fields. In particular, the 
nonlinearity appearing in the two-photon state is attracting much attention due to its possible application to quantum 
information devices and to recent progresses in the photon-pair manipulation technique 

In order to discuss the dynamics of two-photon state theoretically, quantum-mechanical treatment of the photon 
field is indispensable, because the nonlinearity appears in the wavefunction of the photon, not in the amplitude Q. 
■ Such an analysis was pioneered in a case where the nonlinear material is a one-dimensional atom, which is obtained 
as the lossless and weak-coupling limit of the atom-cavity system ||. However, considering that the photon field is 
more magnified for better (higher Q-value) cavities, the strong-coupling cases also seem promising in yielding large 
O |' nonlinearity The optimum cavity conditions for the two-photon nonlinearity have not been sufficiently discussed, 
which should be studied with the method applicable to any coupling regime. In this study, we derive the analytic 
expression of the two-photon propagator in general atom-cavity system, where both weak and strong coupling regimes 
are covered and the loss from the cavity due to radiative decay into non-cavity modes is taken into account. Using this 
propagator, we discuss the two-photon nonlinearity, and clarify the optimum condition on the photonic wavepacket 
for achieving large nonlinearity depending on the cavity conditions. 

The composition of this paper is as follows. In Sec. [HI the theoretical model of the atom-cavity system is introduced, 
and the meanings of the parameters are described. In Sec. IIIII we define the input and output states of photons. In 
Sec. IIVI the measure of the nonlinearity appearing in the output state is defined. In Sec. we show the form of the 
input wavefunction and a scaling law for this atom-cavity system. In Sec. lVIl we numerically evaluate the nonlinearity 
for the lossless case, and clarify the optimum condition for inducing large nonlinearity. In Sec. IV11I the effect of the 
loss is discussed. The analytic expressions for the one- and two-photon propagators are shown in Appendix. 

II. SYSTEM 

As the nonlinear optical system, we here investigate a single two-level system (hereafter called an "atom") embedded 
in a one-sided cavity ^(j ■ The system is schematically illustrated in Fig. 2] The atom is coupled, besides the cavity 
mode, to the photon field in the lateral direction. The cavity mode is coupled through the right mirror of the cavity 
to the external photon field, which is labeled one-dimensionally by r. Although the external field actually extends 
only in the r > region and the incoming and outgoing photons are flyin g in the opposite direction, we treat the 



incoming photons to propagate in the r < region in the positive direction 
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FIG. 1: Schematic view of the atom-cavity system. The right mirror of the cavity is weakly transmissive, through which the 
cavity mode is coupled to the external photon field, g, 7, k represents the atom-cavity coupling, the radiative decay rate of the 
atom into lateral photon modes, and the decay rate of the cavity mode, respectively. 



The Hamiltonian of the system is given, putting K = c = 1, by 

H = u> a a + a- + Jdn^dj l d tJ ,+ J dpi (j^a+d^ + h.c.) 
+ u) c c^c + J dk kb\.bk + J dk (Kkc'bk + h.c.) 

+ g(a+c + Ja-), (1) 

where <r_, c, bk, and d^ are the annihilation operators for the atom, cavity mode, external photon mode, and lateral 
photon mode, respectively. uj a and lo c represent the frequencies of the atomic transition and the cavity mode, and g 
represents the atom-cavity coupling. Regarding 7^ and Kk, we use the flat-band assumption, i.e., 7^ = (7/2^) x / 2 and 
Kk = (K/27T) 1 / 2 , by which the damping rates of the atom and the cavity mode are given by 7 and k. 
The complex frequencies of the atom and cavity are defined by 

u} a =uj a - 17/2, (2) 
uj c = lo c — in/2. (3) 

Using these frequencies, the complex eigenfrequencies d>i 2 of the atom-cavity system are given by 

(uj - uji)(uj - uj 2 ) = (w - u a ){uj - Q c ) - g 2 . (4) 

We remark here that 7 is related to the dissipation of this atom-cavity system. When 7 = 0, the inputted photons 
are always reflected back into the output port, i.e., the atom-cavity system is lossless. Contrarily, when 7^0, some 
of the inputted photons are lost as the spontaneous emission into the lateral direction, i.e., the atom-cavity system 
is lossy. We also remark that, in the limit of 7 — * and (g, k) — > 00 keeping T = Ag 2 /K constant, this atom-cavity 
system is reduced to the one-dimensional atom, where the atom is directly coupled to the external field at r — with 
the coupling constant T 0,0- 

The real-space operator of the external field, b r , is given by the Fourier transform of 



b r = {2-k)- 1 ' 2 j dke ikr b k . (5) 
We note again that the negative (positive) r represents the incoming (outgoing) fields. 



III. INPUT AND OUTPUT STATES 



Our main concern in this study is how the initial one- or two-photon wavepackets are transformed after the inter- 
action with the atom-cavity system. In the initial state, the atom-cavity system is in the ground state and one or two 
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photons exist in the input port, i.e., 

I*^> - / dnMr)&t|0>, (6) 
l*L 2) ) = 2- 1 / 2 j dr 1 dr 2 ^„(r 1 ,r 2 )6t i 6tj ). (7) 

The one-photon wavefunction tpi n (r) is normalized as / dr\^p- ln {r)\ 2 = 1, and ipin(r) = in r > 0. Similarly, Vin(^i, ^2) 
satisfies J <ir 1 c?r2 1 "0in ( r 1 5 r 2)| = 1 an d ^inC^i) ^2) = in ri,r 2 > 0, and has the symmetry ip'm(r 2 , — V'in(Pi) ^2)- 

Sufficiently after the photons have interacted with this atom-cavity system, the excitations in the atom and the 
cavity mode completely escape to the external modes or to the lateral modes Then, the states of the system 
are written as 

l*&> = / dWw(r)fot| } + ..., (8) 
|*£!> - 2-V2 J dr 1 dr 2 ^ out (r u r 2 )bt i bt 2 \0) + ---, (9) 

where the dots imply the terms containing the lateral mode excitations. While the number of the photons are kept 
in the output wavepacket (J dr\tp out (r)\ 2 = 1 and J dridr^Vw^i, ?*2)| 2 = 1) in case of 7 = 0, the output state is 
attenuated in general cases of 7 ^ 0. 

The input and the output wavefunctions are related by the propagator G as follows; 

Vout(r) = J dr'G(r;r')i> in (r'), (10) 

Vw(ri,r 2 ) = J dr' 1 dr 2 G(ri,r 2 ;r' 1 ,r 2 )ip in (r' 1 ,r 2 ). (11) 

The one- and two-photon propagators are analytically obtainable for this atom-cavity system, which are shown in the 
appendix. The two-photon propagator has the symmetry of G(r 2 , r\; r' 2 , r[) — G(n, r 2 ; r[ , r 2 ), which guarantees the 
symmetry in the output wavefunction, Vw( r 2j r i) = Vw^i, r 2 ). 

IV. MEASURE OF NONLINEARITY 

When two photons are inputted into this atom-cavity system, the input wavefunction fAm^i, ^2) is finally trans- 
formed to the output wavefunction ipout(ri 7 r 2 ). In order to quantify the nonlinear effect in this process, we compare 
Vw^i) ?*2) with the linear output wavefunction ^ ut (ri,r 2 ), which is defined by 

Ct(ri,r 2 ) = |drWr 2 G(r i; ri)G(r 2 ;r 2 )Vi„(ri,r 2 ), (12) 

where G(r; r') is the one-photon propagator. This linear output is obtained when the atom in the cavity is replaced 
by a harmonic oscillator with the same transition frequency, i.e., when the nonlinearity of the system is completely 
removed. 

We here define the following quantity (3 by 

= /^idr 2 (V'out)*V'out ^ 

,J J c?ri(2r 2 |Vw| 2 \J J c?ridr 2 |i/>out| 2 

(3 always lies in the unit circle (|/3| < 1) due to the Schwartz's inequality, and (3 = 1 holds when the response of the 
system is completely linear (Vwt = V'out)- Thus, the nonlinear effect is reflected in the deviation of (3 from the unity, 
and \j3 — 1 may be regarded as a measure of the nonlinear effect. We also note that, when the atom-cavity system 
is lossless (7 = 0), the norms of Vw and Vout i n tne denominator are always unity, and (3 is simply reduced to the 
overlap integral between Vw and V'out- 
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V. INPUT WAVEFUNCTIONS AND SCALING LAW 

Now we embody the above formalism for specific forms of the input states, and clarify the conditions to obtain 
large nonlinearity between two photons. In this study, we focus on the case where the two input photons have the 
same wavefunction of the Gaussian form, i.e., 

ipin(ri,r 2 ) = i>m{ri— a) ip- m (r2 — a), (14) 
ip in (r) = (2/ 7 rd 2 ) 1 / 4 exp(-r 2 /d 2 + i gr). (15) 

The input wavefunction is characterized by two parameters, q (the central frequency) and d (the coherent length of 
the photon). a(< 0) in Eq. (|14J) represents the initial position of the photons, which is an irrelevant parameter. 

The atom-cavity system is characterized by g (the coupling between the atom and the cavity mode), uj a — lo c 
(the frequency mismatch between the atom and the cavity mode), 7 (the damping rate of the atom into the lateral 
modes) , and n (the damping rate of the cavity mode) , Adding the input-photon parameters q and d, the whole system 
is specified by the following set of parameters, {g,0J a — w c ,7,K,g,d). However, the number of the parameters may 
be decreased with the help of the scaling law that the system with the parameters (ag, a(w — u> c ), cry, an, aq, d/a) 
is equivalent to the one with {g,oj a — w c , 7, k, q, d), where a is a positive constant. Using this law, the system is 
specified by the following set of dimensionless parameters, ((w a — uj c )/g, 7/(7, k/<?, q/g, gd). uj a — lu c is fixed at zero in 
the following numerical examples, and uj a {—uj c ) is chosen as the origin of frequency. 

VI. NUMERICAL RESULTS FOR LOSSLESS CASES 

In this section, we show the numerical results for 7 = cases, where the atom-cavity system is lossless and the 
number of photons is preserved in the output wavepacket. The atom-cavity system is then characterized only by the 
ratio n/g, as has been discussed in Sec. E] Eq. 0} indicates that the Rabi splitting of the eigcnfrcqucncy of the 
atom-cavity system takes place when n/g < 4. Thus, in our definition of the parameters, the strong (weak) coupling 
regime is specified by n/g < 4 {n/g > 4). 

A. Weak coupling regime 

First, we discuss the weak coupling regime. In Fig. the nonlinearity \j3 — 1| is plotted for the cases of resonant 
input (q — 0) as a function of k/ g and g 2 d/n. It has been confirmed that the nonlinearity appears most strongly at 
q = in the weak coupling regime and the nonlinearity gets weaker as |g| is increased. 

Fig- [21 indicates that, roughly speaking, the nonlinearity \/3 — 1| is dependent solely on g 2 d/n in the weak coupling 
regime (n/g > 4). This is because, in the weak coupling regime, the atom-cavity system may be regarded as a "one- 
dimensional atom" , where the atom is coupled directly to a one-dimensional photon field with a coupling constant 
4<? 2 /k. It is observed that the nonlinearity is maximized when g 2 d/n ~ 0.5. For example, for g =120MHz and 
k=900MHz , the optimum pulse length d is about 9m. The qualitative explanation for this optimum condition will 
be given in Sec. IVI CI 

B. Strong coupling regime 

Next, we discuss the strong coupling regime. In Fig. the nonlinearity \0 — 1| is plotted for the fixed «/<?(=(). 5) 
as a function of nd and q/g. Contrarily to the weak coupling case, the nonlinearity is weak for the resonant input 
(q = 0). This is because the q = photons are no more resonant to the cavity mode due to the Rabi splitting. Instead, 
the nonlinearity appears strongly when the input photons are tuned to the Rabi-splitted resonant frequency, q ~ ±g. 
The maximum value of \[3 — 1| is about 1.5, which is almost the same as the weak coupling case. The nonlinearity is 
optimized for nd ~ 4, the reason of which will be discussed in Sec. IVI CI 

C. Optimum pulse for inducing large nonlinearity 

In the preceding subsections, we have clarified the optimum frequency q and the length d of the photon pulse for 
inducing large nonlinearity, in both weak and strong coupling regimes. Here we explain these optimum conditions 
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FIG. 3: Dependence of \/3 — 1| on nd, where n/g is fixed at 0.5 (strong coupling regime), q/g = 0.8 (broken line), q/g — 0.9 
(solid line), q/g = 1 (dotted line), and q/g = (thin broken line). 



from a unified viewpoint. To this end, we focus on the following wavefunction y(r), which is defined by 

e mt a + \0) = J dnp{r)b\.\Q) (16) 

for large t(> 0). The meaning of <p(r) becomes clear by multiplying e~ lHt from the left; if a single photon pulse (p(r) 
is prepared as an input at time to, the photon will completely be absorbed by the atom at time to + 1. Therefore, it 
is expected that, when the input pulse V , in(?") resembles f(r) in shape, the two input photons try to occupy the atom 
simultaneously and the nonlinearity appears strongly. 
<p(r) has the following form; 

,s \ i^l( e ^T (»-+*)- e l -2 (r+t)) (-t< r <0) 
(otherwise) 

where uj\ and UJ2 are the complex eigenfrcquencies of the atom-cavity system, defined in Eq. In the weak coupling 
regime, are approximately given by —in/2 and —2ig 2 /n. Noticing that g 2 /n ^ k in this regime, the optimum 
frequency q and pulse length d are given by q ~ and d ~ n/2g 2 , which explains the numerical results in Sec. IVI 
On the other hand, in the strong coupling regime, Cj\^ are approximately given by — in/A±g. Therefore, the optimum 
q and d are roughly estimated at q ~ ±g and d ~ 4/k, which are in agreement with the numerical results in Sec. IVI 51 
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FIG. 4: Norms of ipout{ri, f"2) (bold lines) and V'out('"i, r 2) (thin lines), for lossy cases of j/g = 0.1 (solid lines) and j/g — 0.1 
(broken lines). The atom-cavity system is in the weak coupling regime (n/g — 5) and the resonant photons (q = 0) are used. 
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FIG. 5: Effect of the loss on the nonlinearity |/3 — 1|; 'y/g = (thin dotted line), 'y/g = 0.1 (solid line), and f/g = 0.2 (broken 
line). Other parameters are nj g = 5 and q — 0. 



VII. NUMERICAL RESULTS FOR LOSSY CASES 

In the preceding section, the results for the lossless cases (7 = 0) are presented. Here, we discuss the lossy cases 
(7 / 0), fixing the parameters n/g — 5 (weak coupling regime) and q = for example. 

The apparent result of the presence of the loss is attenuation of the photon pulses. In Fig.0J we have plotted the 
norm of the two-photon output wavefunction ■0out( r i; r 2)i i-e. , the probability to find two photons in the output. The 
norm of the linear output wavefunction ^ ut (ri, r 2 ) is also plotted in the same figure. As expected, the photon pulse 
is attenuated more significantly when the loss parameter 7 is larger, for both i/> ou t and V'out (compare solid and broken 
lines). It is also observed that f/'out * s m ore attenuated than ip ou t- This is understood by the facts that, in the linear 
case, the photons are more likely to be absorbed by the atom due to the absence of the saturation effect, and that 
the loss of photons may take place only while the photons are occupying the atom, as schematically shown in Fig. ^ 

Fig- El plots the nonlinear parameter \(3 — 1| defined in Eq. 113|) . when the system has finite loss 7. It is observed 
that |/3 — 1| is slightly decreased (increased) for g 2 d/ n < 0.5 (g 2 d/n > 0.5) region but, overall, that no qualitative 
changes are introduced by the presence of the loss. Thus, although the probability to find two photons in the output 
is largely decreased as shown in Fig. ^ the nonlinear characteristics of the output state are almost unchanged from 
the lossless case if the output pulse contains two photons. 
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VIII. SUMMARY 

We have theoretically investigated a situation where two photons are simultaneously inputted into a nonlinear 
optical system, and examined the nonlinear effect in the output wavepacket. As a model nonlinear optical system, we 
employed a two-level system (atom) embedded in a cavity, which is illustrated in Fig. ^ The one- and two-photon 
propagators are obtained in analytic forms, which are presented in Appendix. The overlap integral of the linear 
and nonlinear two-photon output wavepackets, Eq. (|13l) . is used to measure the nonlinear effect. This quantity is 
numerically evaluated both in the weak coupling regime (Fig. [5J and in the strong coupling regime (Fig. EJ, and the 
conditions for inducing large nonlincarity is revealed. These conditions arc explained by the optimum pulse shape, 
Eq. (|17l) . This finding suggests that the pulse shape control is more essential in optimizing the two-photon nonlinearity 
rather than the Q-value control of the cavity system. We also considered the 7^0 cases, where the atom-cavity 
system is lossy. As shown in Fig. 0J the probability to find two photons in the output is largely decreased by the effect 
of the loss, but the nonlinear characteristics of the output two-photon output state are almost unchanged from the 
lossless case. 
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APPENDIX A: PROPAGATOR 

In this section, we present the analytic forms of the one- and two-photon propagators G(r; r 1 ) and G(r%, r 2 ; r^r^). 
Using the completeness relation in the one-photon space, 1 = J drbl\0)(0\b r , the one-photon propagator is given by 

G(r;r') = {0\b r e- mt bi,\0), (Al) 

where r > > r', and t is the time difference between the output and the input. This propagators are derivable 
by standard application of the Green function method [T^ | . Throughout this section, we use the coordinate system 
moving at the light speed; transformation to the static coordinate system is done by replacing r (coordinates without 
primes) by r — t, or r' (coordinates with primes) by r + 1. In this coordinate system, the one-photon propagator is 
identified as 



G(r, r') = G(r - r') = G (r - r') + G 2 (r - r'), (A2) 

Go(r-r') = S(r -r') ~ K0(r' ~r)e iCj ^ r - r '\ (A3) 

G 2 {r - r') = K6(r' - r) ( j^-r>) _ ^ ~ "a j Ql (r-r') _ "l-"c e i£ 2 (r-r'A (A4) 

V C02 — 0J\ UJ\— UJ 2 / 



where uj a (complex frequency of the atom), lu c (complex frequency of the cavity), and d>i 2 (complex eigenfrequencies 
of the atom-cavity system) are defined in Eqs. J3J, (J3J), and (J2J. 

Next, we proceed to the two-photon case. Using the completeness relation in the two-photon space, 1 = 
2 _1 J dridr2bl 1 bl 2 \0)(0\b ri b r2 , the two-photon propagator is identified as 2~ 1 (0|& ri 6 r2 e~ lW *^, b' r , |0). This quantity 

is composed of two kinds of terms; in the first (second) kind, the photons initially located at r[ and r' 2 are scattered 
to n and r 2 (r 2 and n), respectively. With the help of ipi n (r[ , r' 2 ) = ipi a (r' 2 , r'{), it is shown that both of them yield 
the same output wavefunction V'out('"i,'"2) after integration in Eq. Qlip. We can therefore safely regard only the first 
kind of terms as the two-photon propagator. It is given by 

G(r 1 ,r 2 ,r' 1 ,r 2 )=G(r 1 -r' 1 )G(r 2 -r 2 )-G 2 (r 1 -r' 1 )G 2 (r 2 -r' 2 )+ £ G^n,^,^), (A5) 

J=4,6,8 

where Gj (j — 4, 6, 8) are defined by 

G A {r u r 2 ,r' 1 ,r' 2 ) = - [J 4 (n - r 2 , r[ - r' 2 , r 2 - r[) + h{r 2 - n, r 2 - r[, r x - r' 2 )] , (A6) 
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where 



G 6 (ri,r 3 ,r' u r' 2 ) = 3 [h(r 2 - n, r[ - r' 2 ,n - r[) + h(r x - r 2) r' 2 - r[,r 2 - r' 2 )] , (A7) 
Gs^,^,^,^) = ~^ 3 K [J 8 (n - r 2 , r\ - r' 2 , r 2 - r[) + J 8 (r 2 - n, r' 2 - r[,n - r' 2 )\ , (A8) 

h(x,y,z) = dkdqdu—. J{ "\ q ' K) , (A9) 

7 (w — A; — g + id) 

I 6 (x,y,z) = / dkdqduj- — — — — — , (A10) 

T< s /■,,,, e^+^y+^J{LO,q,k) 

h{x,y,z) = / dkdqduj- — — — — — , (All) 

J {u - k - lu c )(u) - q - u c )(uj - Vq)(w - - v 2 ) 



J (w, q, k) = — — — — — — — (A12) 

(k ~ bJ c ){us — k — u>i){u> — k — uj 2 )(q — oj c )(uj — q — u>i){u> — q — uj 2 ) 

and ^0,1,2 are defined by (to — vo){uj — Vx){u> — v 2 ) = (lo — uj a — uj c )[{u — 2u! c )(uj — ui a — uj c ) — 2g 2 ]. The triple integrals 
in the definitions of ^4,6, 8 are carried out analytically. It can be confirmed that the nonlinear part of the Green 
function, G NL = -G 2 {r\ - r' 1 )G 2 {r 2 - r 2 ) + J2 J= 4,6,s Gj(ri,r 2 ,r[,r 2 ), is nonzero only when max(ri, r 2 ) < min(ri, r' 2 ) 
is satisfied, i.e., the nonlinearity appears only after the two photons have arrived at the atom. The one- and two- 
photon propagators reduce to those for the one-dimensional atom in the limit of 7 — > and (g, k) — » 00 keeping 
T = 4g 2 /« constant @. 
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